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ON  MINIMAX  ROBUST  DATA  FUSION 
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and  System  Research  Center 
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College  Paris:,  MD  20742 


ABSTRACT 

In  this  paper,  minimax  robust  data  fusion  schemes  based  on  discrete-time  observations 
with  statistical  uncertainty  are  considered.  The  observations  are  assumed  to  be  i.i.d  and  the 
decisions  of  all  sensors  independent  when  conditioned  on  the  either  of  two  hypotheses..  The 
statistics  of  the  observations  are  only  known  to  belong  to  uncertainty  classes  determined  by  2- 
altemating  Choquet  capacities.  Both  cases  of  fixed-sample-size  (block)  data  fusion  and  sequen¬ 
tial  data  fusion  are  examined.  For  specific  performance  measures,  three  robust  fusion  rules: 
suboptimal,  optimal  and  asymptotically  optimal  --as  the  number  of  sensors  increases-are 
derived  for  the  block  data  fusion  case,  and  an  asymptotically  robust  fusion  rule  is  derived  for 
the  sequential  data  fusion  case;  these  fusion  rules  are  optimal  in  the  class  of  rules  employing 
likelihood  ratio  tests.  In  all  situations  the  robust  fusion  rule  makes  use  of  likelihood  ratios  and 
thresholds  which  depend  on  the  least-favorable  probability  distributions  in  the  uncertainty  class. 
In  the  limit  of  a  large  number  of  sensors,  it  is  shown  that  the  same  threshold  can  be  used  by 
all  sensors,  which  in  turn  simplifies  the  overall  computation. 


This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under  contract  N00014-86-K-0013  and  in  part  by  the 
Systems  Research  Center  at  the  University  of  Maryland,  College  Park,  through  the  National  Science  Foundation's  Engineering 
Research  Centers  Program:  NSF  CDR  8803012. 


I.  INTRODUCTION  AND  PROBLEM  FORMULATION 


Recently,  distributed  detection  with  or  without  a  fusion  center  has  attracted  consid¬ 
erable  attention([l]-[3]).  In  [1]  it  is  shown  that  the  optimal  decision  statistic  for  fixed- 
sample-size  (block)  detection  is  a  likelihood  ratio  test  with  dependent  thresholds  which 
satisfy  coupled  equations.  The  computation  of  these  thresholds  is  a  very  difficult  task 
and  becomes  even  more  difficult  when  sequential  test  is  used  ([2]).  Further  in  [l]-[3]  the 
statistics  of  the  observations  are  assumed  to  be  known  a  priori,  which  can  be  unrealistic 
in  practical  situations.  Frequently  only  imperfect  knowledge  is  available  in  real  situations. 
In  [4]  minimax  robust  distributed  detection  schemes  without  fusion  were  developed  for  a 
general  uncertainty  class,  the  class  of  2-alternating  Choquet  capacities.  In  this  paper  min¬ 
imax  robust  data  fusion  schemes  are  derived  for  the  binary  discrimination  problem.  The 
uncertainty  class  considered  here  is  the  same  as  that  of  [4],  it  contains  many  uncertainty 
models  ([5]- [7])  and  is  described  by 

V  =  {P|P(S)  <  v(S),  VS  €  B,  P(ft)  =  r(ft)}  (1) 

where  ft  is  a  sample  space,  B  the  associated  cr-field  and  v  a  2-alternating  capacity  [5].  The 
properties  relevant  to  our  paper  of  this  uncertainty  model  are  stated  in  the  following. 
Lemma  1:  Suppose  and  Vi  are  two  alternating  capacities  on  (ft,P)  and  Vq  and  V\ 
are  the  uncertainty  classes  determined  by  them  as  in  (1).  Then  there  exists  a  Lebesgue- 
measurable  function  it  :  ft  — ♦  [0,  oo]  such  that 

9vq({tcv  >  9})  +  »!({*„  <  9})  <  9v0(S)  +  t>! ( Sc)  (2) 

for  all  S  G  B  and  all  9  >  0.  Furthermore  there  exist  measures  (Pq,  A)  in  V,  X  V\  such 
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that 


A({jt»»  >  0})  =  Vo({7r„  >  0}) 

(3) 

A({*\.  <  9})  =  Vl({7T„  <  ^}) 

(4) 

i.e.,:  t„  is  stochastically  largest  over  Vo  under  Pq  and  stochastically  smallest  over  V\ 
under  Pi  .The  function  irv  is  a  version  of  dPi/dP0  and  unique  a.e.  [Po  + A]-  The  measures 
(Po,  A)  are  termed  the  least-favorable  measures  over  Vo  x  V\. 

The  system  models  used  here  are  similar  to  those  of  [2]  and  [9].  Each  sensor  makes  a 
decision  based  on  the  observations  i /,(<)  and  sends  a  binary  message  U{  (t  =  1,2,  •  •  •,  I() 
which  represents  its  decision  ;  then  the  fusion  center  makes  the  final  decision  according  to 
a  specific  fusion  rule.  Hereafter  the  following  assumption  is  made. 

Assumption  1(A1):  (i)  The  statistics  of  the  observations  conditioned  on  either  hypothe¬ 
sis,  Pi(yj(t))  (i  =  0, 1;  j  —  1, 2,  •  •  • ,  K),  are  the  same  for  each  sensor  j  and  all  time  instants 
t.  (stationarity)  (ii)  p, •(%(<))  are  independent  across  sensors  and  time  instants. 

Both  cases  of  fixed-sample-size  (block)  data  fusion  and  sequential  data  fusion  are 
considered.  Denote  the  likelihood  ratio  and  the  optimal  threshold  of  sensor  i  by  Z,  and 
rji,  respectively,  ( i  =  1,  -  *  • ,  J5f ),  and  by  rjo  the  optimal  threshold  of  the  fusion  center.  In 
the  case  of  block  data  fusion,  when  K  is  finite,  the  performance  index  adopted  here  has 
the  form 

CK{Li,rji 5  *  =  0 

=  c\PF0{Li,rn\i  =  0,---,tf)  +  c2PM0(T,'»7?,-;t  =  0,  -  -  - ,  K)  (5) 

where  Pp0  and  Pm0  are  the  probabilities  of  false  alarm  and  missing  of  the  fusion  center  and 
ci,c2  are  positive  constants  representing  cost.  Further  let  pm(n*  =  0,1)  denote  the  prior 
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probability  of  hypothesis  Hm  then  Pt  =  PqPf0  +  Pi  Pm0  is  the  error  probability.  When 
K  — *•  oo,  since  the  error  probability  will  approach  zero,  the  performance  is  measured  by 


Joo{Li,  Tfi  —  0, 1,  ■ 


• , )  =  lim 

’  K-+  oo 


In Pe(Lj,rji,  i  =  0,---,K) 
K 


(6) 


which  is  the  exponential  rate  of  decrease  of  the  error  probability  and  is  useful  for  asymp¬ 
totic  analysis  ([2]&[8]).  For  sequential  data  fusion,  besides  the  error  probability,  the 
expected  value  of  the  sample  size  of  the  fusion  center  must  also  be  included  as  part  of  the 
performance  measure. 

This  paper  is  organized  as  follows.  In  section  II,  robust  block  data  fusion  algorithms 
are  derived  for  suboptimal  fusion  and  for  optimal  fusion  with  equal  local  thresholds.  In 
particular  as  the  number  of  sensors  K  —*  oo,  the  asymptotically  robust  fusion  rule  is 
derived.  In  section  III,  an  asymptotically  robust  sequential  data  fusion  rule  is  developed. 
Finally, 
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II.  MINIMAX  ROBUST  BLOCK  DATA  FUSION 


II. A  “AND”,  “OR”  and  Optimal  Fusion  Rules 

Two  simple  data  fusion  rules  are  :  (1)  The  fusion  center  decides  H\  if  all  sensors 
transmit  u,  =  1  for  all  i  — termed  “AND”  rule.  (2)  The  fusion  center  decides  H\  if  at 
least  one  sensor  transmits  «,•  =  1  for  1  <  t  <  K  — termed  “OR”  rule.  Notice  that  for 
these  two  suboptimal  fusion  rules,  there  is  no  t?o  in  the  usual  sense  of  thresholds.  These 
two  policies  are  suboptimal  but  easy  to  implememt,  especially  when  the  number  of  sensors 
become  very  large. 

Under  (Al),  the  error  probabilities  of  false  alarm  and  missing  can  be  expressed  as 

K 

PFo  =  P(  m  =  1  ;  t  =  1,  •  •  • ,  K\H0)  =  n  PF  (7) 

i=i 

and 

K 

PMo  =  P(vi  =  <m<i<K\H1)  =  Y,PM‘  (8) 

«=i 

for  the  “AND”  fusion  rule;  and  as 

K 

PFo  =  P(ui  =  1  ;1  <  i  <  K\Hq)  =  '£PF  (9) 

t=i 


Usually  for  the  purpose  of  performance  evaluation  the  distribution  of  the  observations 
is  assumed  to  be  known  a  priori.  However,  what  we  may  have  in  practice  is  some  imperfect 
knowledge  about  the  statistics  of  observations  .  Let  the  statistics  of  data  belong  to  the 
uncertainty  set  defined  by  equation  (1).  Then  by  lemma  1  the  error  probabilities  PF<  and 
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P\ft  of  each  sensor  (1  <  t  <  K)  satisfy  2  [5] 

PFi(  Li,  r,i )  <  PFi{  Li,  fji)  )  <  PF,(  Di ,  m  )  (11) 

and 

Pm,(  Li,  Vi )  <  Lx,  Vi )  $  Pm,(  Di,  Tji )  (12) 

where  Li  =  dPu/Poi  is  the  worst-case  likelihood  ratio  of  the  ith  sensor  (the  measures  Pu 
and  Poi  are  singled  out  by  lemma  1)  and  Di  is  any  other  potential  test  of  that  sensor  with 
the  same  observed  data.  For  finite  summations  ([8]&[9])  or  multiplications  ([7])  over  i , 
the  direction  of  inequalities  still  holds  the  same  as  in  (11)  and  (12).  Thus  the  following 
proposition  holds. 

Proposition '  1:  Suppose  the  performance  index  of  data  fusion  with  “AND”  or  “OR” 
fusion  rules  is  defined  by  equation  (5).  Then  under  (Al) 

Ck(  Li,  fji",  i =  !,•••,  K ) 

<  Ck(  Li,vr,  i  =  1,---,K) 

<  Ck{  Di,r]i;  i  =  (13) 

Remark  1:  The  threshold  r?o  or  r?0  of  the  fusion  center  are  missing  from  (13)  since  the 

“AND”  or  “OR”  fusion  rules  do  not  involve  common  thresholds. 

Next  we  consider  the  optimal  robust  fusion  rule.  Here  the  optimality  is  in  the  sense 

that  the  center  uses  the  Neyman-Pearson  test  oh  the  transmitted  binary  signals  by  the 

sensors.  When  all  the  sensors  are  alike  and  with  the  same  level  of  error  probabilities,  i.e., 

2 a  variable  with  a  hat  means  that  it  is  a  function  of  worst-case  values;  a  probability  with  a  hat  means 
that  the  corresponding  sensor  operates  under  worst-case  conditions. 
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Pp{  =  Pp  and  PM,  =  Pmi  thus  =  r](i  =  1  the  Neyman-Pearson  test  becomes 

k  —  out  —  of  —  K  rule  [3]  under  (Al).  Therefore  the  optimal  fusion  rule  is 

k  <Ho  k *  (14) 

where  k *  is  a  suitable  threshold  related  to  t)q.  In  this  case  the  overall  error  probabilities 
of  false  alarm  and  missing  under  mismatch  [that  is,  k *  is  the  threshold  matched  to  the 
worst-case  ( Pf,Pm ),  but  (Pp,Pm)  characterize  the  operating  conditions]  are 

PFo=  £  (fW(l-PF)(*-‘)  (15) 

«'=ffc*l  \  / 

Pm o=  £  (^J(1-.Pm)‘V*-,)  (16) 

Notice  that  A:*  is  a  function  of  Pp,  Pm  and  K  ;  Pp  and  Pm  are  functions  of  the  number 
of  samples  n  that  each  sensor  collects  for  making  its  decision.  For  proposition  2,  another 
assumption  is  made  (see  Appendix) 

Assumption  2(A2):  Pp  <  \k*]/K  and  Pm  <  1  -  [A:*] /K  for  all  K. 

Then  from  the  above  definition  of  Pp0  and  Pm0,  proposition  2  follows. 

Proposition  2:  Suppose  the  thresholds  of  all  sensors  are  equal,  then  under  (Al)  and 
(A2)  the  following  inequalities  are  true 

Cp{  Li,fj ;  i  =  0,-  •  • ,  K) 

<  Cp{  -£«,  *);  *  =  >  k) 

<  CK{Di,T1-,i  =  0,...,K)  (17) 

Proof :  Taking  derivative  of  (15)  with  respect  to  Pp ,  we  have 

^r  =  £  (^P^ii-PF^-^ii-KPp) 

F  *=r*-i v 
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If  (A2)  is  satisfied,  then  Pf  <  Pf  <  \k*']/K  as  well  and  subsequently  dPpQ /dPp  >  0 
;thus  Pp0  is  an  increasing  function  of  Pp -  By  a  similar  argument,  it  can  be  shown  that 
dP\tol dPiA  >  0,  thus  Pm0  is  also  an  increasing  function  of  Pm.  Hence  Cp  will  be  an 
increasing  function  of  Pf  and  Pm  under  (Al)  and  (A2).  Thus  by  (11)  and  (12),  (17) 
holds. 

Remark  2:  The  assumption  (A2)  is  nonrestrictive  for  practical  cases  (refer  to  Appendix). 
If  it  is  not  satisfied,  data  fusion  is  of  no  use  because  Pp0  and  Pm0  do  not  offer  any 
improvement  over  Pp  and  Pm- 

The  above  proposition  only  covers  a  special  case  [9];  in  general  the  computation  com¬ 
plexity  of  the  optimal  thresholds  of  the  sensors  is  N-P  complete.  Thus  our  interest  in  the 
limiting  case  ,  K  — ►  oo,  is  justified. 

II.B  Asymptotic  Behavior 

When  the  number  of  sensors  increases  to  infinity,  the  probability  of  error  decreases 
to  zero.  Therefore  the  performance  index  adopted  in  asymptotic  analysis  should  be  the 
exponential  rate  of  decrease  of  the  error  probability  (6),  where  P£  =  poPp0  +  P\Pm0- 
As  stated  in  theorem  1  of  [2],  when  K  — ►  oo  under  the  assumption  1  of  [2],  the  same 
threshold  r]  could  be  used  in  each  sensor  for  optimal  data  fusion  (to  minimize  In  Pe/I()). 
Hence  J0 Q(  Dj,T){;i  =  0,  l,---)  =  /<„(  D{,v,  i  =  0,  l,*-*)  and  the  error  probabilities  Pp0 
and  Pm0  satisfy  (15)  and  (16).  Then  proposition  3  follows. 

Proposition  3:  Under  (Al)  and  (A2),  as  K  — ►  oo,  the  performance  index  defined  by  (6) 
satisfies 


Joo{  Vo,  V,  *'  —  6, 1,  •  •  ■) 
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<  Joo(Lit  fK>,fi;i  =  0,1,-  ••) 

—  «/oo(  »?0,  *?!  *  =  0,  1,  ■  •  •) 


(18) 


Proof:  Based  on  the  above  discussion,  a  similar  proof  to  that  of  proposition  2  can  be 
developed  by  replacing  Ck  with  Jqo- 
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III.  ASYMPTOTICALLY  ROBUST  SEQUEN¬ 
TIAL  DATA  FUSION 

In  this  section  the  SPRT  (sequential  probability  ratio  test)  is  adopted  as  the  decision 
policy  of  the  fusion  center  whereas  the  policies  of  all  the  sensors  remain  the  same  ,  i.e., 
likelihood  ratio  tests.  The  data  used  for  fusion  are  the  binary  signals  transmitted  by  the 
sensors. 

The  scheme  of  sequential  data  fusion  is  used  to  minimize  the  following  cost  function  : 
the  exponent  of  the  sum  of  the  error  probability  and  the  average  sample  sizes  under  both 
hypotheses.  This  kind  of  fusion  rule  has  not  yet  been  shown  to  be  the  optimal  one,  but  for 
robustness  in  the  class  of  likelihood  ratio  tests, it  is  a  reasonable  choice  for  the  sequential 
fusion  rule.  Thus,  in  this  case  ,  the  expected  value  of  the  sample  size  No  of  the  fusion 
center  (that  is,  the  number  of  times  the  fusion  center  collects  the  decisions  of  the  sensors 
before  reaching  its  decision)  enters  into  the  performance  measure. 

Again  let  the  statistics  of  observations  conditioned  on  either  hypothesis  belong  to  the 
uncertainty  model  defined  by  (1).  When  K  — ►  oo,  under  suitable  conditions,  it  can  be 
shown  that  the  worst-case  pair  of  statistics  of  observations  is  still  the  one  singled  out  by 
lemma  1,  if  the  performance  index  for  robustness  has  the  following  form 

•foo (L{,  TJi',  i  —  0, 1,  •  •  •) 

l 

=  lim/e_co  In  {  y  )  Pm Em {iVo|Lj,  77,;  i  =  0,  •  •  • ,  K }  1 

m= 0 

+Pe(Li,  i)i',  i  =  0,  •  •  • ,  K)}  (19) 

where  rj0  =  (A,  B ),  (A  <  1  <  B)  are  the  optimal  thresholds  of  the  SPRT  used  by  the  fusion 
center  and  77,, i  =  1,  •••,/£'  are  the  optimal  thresholds  used  by  the  sensors.  As  K  — ►  00, 
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the  expected  sample  size  approaches  1  and  the  error  probability  approaches  0.  Thus  (19) 
is  the  exponent  of  the  sum  of  these  two  terms  (  i.e.,  the  sum  is  of  the  form  exp(JooK)) 
and  is  an  extension  of  the  performance  measure  in  [2].  In  the  following,  we  abbreviate 
=  0,---yK)  as  ^  where  2  -'Oh*  V2,-",Vk)  and  =  0  •  *  • ,  A')  as  77  for 

those  expressions  inside  the  conditioning  in  the  expected  sample  size  and  inside  the  error 
probability.  Define  for  0  <  s  <  1 

Mi(7?i,si)  =  ln^/to(01-47.x(ir< 

1=0 

where  is  the  probability  mass  function  under  hypothesis  Hm(m  =  0, 1)  of  the  binary 

message  l  G  {0, 1}  that  the  t — th  sensor  sends  to  the  fusion  center;  /jm(.)  is  stationary, 
independent  of  time.  Let  (77,3)  minimize  /x,(77,-,Si)  for  i  =  1,2,-  •  -(i.e.,  71(77, 3)  <  71;  (77,-,  3,) 
for  all  t)  and  r  =  e**.  Before  we  proceed,  we  state  the  assumption  1  of  [2]. 

Assumption  3(A3):  (i)  174(77,  s)|  <  00,  Vr7,s.  (ii)  There  exists  a  constant  R  such  that 
1/4(77,  s)|  ^  -Ri  V(t7,s)  G  [0,oo]  X  [0,1],  where  a  dot  stands  for  differentiation  with  respect 
to  s. 

'  Notice  that  the  function  /z  has  the  same  form  as  that  in  [2]  and  [10].  In  fact,  we  can 
obtain  a  similar  result  to  theorem  1  of  [2]  for  (19).  First  we  derive  the  upper  and  the  lower 
bounds  on  the  expected  sample  size. 

Lemma  2:  Suppose  (Al)  and  (A3)  are  satisfied  ,  then  for  m  —  0, 1 

+  <2°> 

£m[-Wol2)  >  1  +  ci  exp[A'/z(77,  s)  —  (2KR)1/2}  (21) 

where  do  =  a,di  =  —b  and  ci  is  a  constant  defined  in  the  course  of  the  proof. 
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Proof:  Denote  by  u,j  the  message  send  by  sensor  i  at  time  instant  j  and  by  /;m  its 
(stationary)  distribution  under  hypothesis  Hm{m  =  0,1).  Further  let 


Z,(n)  =  J> 
i= i 


fto(Uij) 


be  the  log-likelihood  ratio  based  on  n  decisions  of  the  i-th  sensor  and  define  by  g(K,n)  = 
K 

y  Zi(n )  the  overall  sufficient  statistic  used  by  the  fusion  center  .Then  the  left-hand  side 
i= 1 

of  (20)  becomes  (a  =  In  A,  b  —  In  B) 


Em[No\v}  =  yPm(No>n ) 

n=0 

oo 

=  1  +  y  Pm(N0  >  n) 

n= 1 

oo  K  * 

=  1  +  ypm(a<yzt(n)<b) 

n= 1  *=1 

OO 

=  1  +  y  Pm(a  <  g(K,n)  <  b) 

n= 1 

Only  the  proof  of  the  case  m  =  0  will  be  given,  since  a  similar  proof  can  be  obtained 
for  m  —  1.  Suppose  the  same  threshold  77  is  used  in  all  sensors,  then 

OO 

y  P0(a  <  g(K,n)  <  b) 

n= 1 

00 

<  y  Po(a  <  g(K,n)) 

n= 1 

n= 1 

=  Ec'”r(vfn 

n=l 

_  e~tar(r},s)K 
1  ~  r{Vi  S)K 

where  the  second  inequality  is  the  Chernoff  bound,  and  the  equality  is  from  the  sum  of  a 
geometric  series  (r(r?,  s)  <  1).  Hence  we  have  (20). 
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On  the  other  hand,  for  the  proof  of  (21),  suppose  different  thresholds  are  used  by 
the  sensors,  then 

> 


where  u.j  are  i.i.d  for  all  j  (across  time).  Let  Uk  =  (a  <  g(K ,  1)  <  6).  Then  by  (3.20)& 
(3.40)  of  [10]  and  under  assumption  1  of  [2] 

Po(UK) 

K  K 

>  c\  exp[^/x,(jjj,Sj)  -  (2^/i1(77,-,si))1/2] 

<=i  ■  i=i 

>  c-iexp[Kfi(T],s)  -  (2KR)1/2} 

where  Qs  and  y,  have  been  defined  in  [10]  and 
0  <  ci  =  YluKny.  Q»(Y)  <  °°.  Thus  we  have  (21). 

Next,  let  us  consider  the  bounds  on  the  error  probability. 

Lemma  3:  Suppose  (Al)  and  (A3)  are  satisfied, 

P'(V)  <  e^(TM)  (22) 

Pe(v)  >  c2exv{Kn(ri,s)E[No]  +  (A>(t7,s))1/2£[JV01/2]}  (23) 

where  E[N0]  =  Em=o  P™E  m  [Ao]  and  c2  is  a  constant  defined  in  the  course  of  the  proof. 
Proof:  Because  A  <  1  <  B, 

P0{L0  >  B\N0  =  n)  <  P0(L0  >  1|A0  =  n); 

P^Lo  <  A[Ao  =  n)  <  Pl(L0  <  1|A0  =  n) 


E P0(a<g(K,n)<b ) 

n=l 

PQ(a  <  g(KJi)  <  b ) 


K 


/«i(u.i) 
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Suppose  the  same  threshold  rj  is  used  for  all  sensors,  then  by  section  III  of  [10],  we  have 


=  f^Pe(v\N0  =  n)P(N0  =  n) 

n= 1 

<  f)  P{N0  =  n)[P0(io  >  1\N0  =  «)  +  Pi(T0  <  l|#o  =  «)] 

n=l 

oo 

<  £  P(No  =  n)eKn^1)'a^ 

n—\ 

<E[eh 


since  p(ri,s)  <  0  and  No  >  1.  Then  we  show  (23).  Denote  by  rj  the  optimal  threshold  for 
block  fusion  under  known  No  =  n,  then 


PoPo(Lq  >  B\N0  =  n)  +  P\P\(Lo  <  AlA^o  =  n) 
=  pqPo(L0  >  v4|iV0  =  n)  +  piPi(Li  <  v4|7V0  =  n) 
>  PoPo(Lo  >  i\N0  =  n)  +  pi  Pi  (To  <  v\N0  =  n) 

According  to  III  of  [10],  we  have 
Pe(v) 

OO 

=  EboPoCzzliVo  =  n)Po(No  =  n) 

n=l 

+PiPi(2l^o  =  n)P\(No  =  n)] 

OO 

>  ^[PoPo(T0  >  rj\N0  =  n)P0(jV0  =  n) 

n=l 

+PiPi(To  <  *i\No  =  n)Pi(No  =  n)] 

oo  1  K 

>  Y2  J2  {PmPm(N 0  =  »)«XpE /*«(»?.•,  Si) 

n=l  m=0  t=l 
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-(2  £>(^*))1/3]  L 

•=1  ymny, 

i  K 

>  Y  c2PmEm{exp[KN0n(q,s)  -  (2 ^  ^t))1^]} 

TO=0  1=1 

where  ^(t  =  0, 1)  were  defined  in  [lOj.  Let 

c2  =  mtn{exp{i;Eyony4gj(y)]},exp{JEEysny>Qa(y)}} 

Using  (A3)  and  the  convexity  of  exponential  function,  we  have 

P'(V)  >  c2exT>{E[KN0fi(r},s)-(2KN0R)1'2}} 

=  c2  expiKfiir,,  s)E[N0]  -  (2 KR)1/2E[N01/2}} 

which  is  (24).  According  to  lemmas  2  and  3,  it  can  be  shown  the  performance  (measured 
by  Joo)  wdl  not  be  hurt  if  the  same  threshold  is  used  by  all  sensors,  and  this  value  could 
be  computed  via  a  simple  method  (optimize  m  with  respect  to  17,  s).  Thus  the  following 
proposition  holds. 

Proposition  4:  Under  (Al)  and  (A3), 

i&o  It  ln{-E'tAr°l2]  “  1  +  p=(l?)} 

=  Urn  ±\n{E[N0\v]-l  +  Pe(v)}  (24) 

K— ►  00  Jti 

Proof:  This  is  a  direct  result  from  (20)-(23)  of  the  above  two  lemmas.  By  using  the  lower 
bounds  in  (21)  and  (23)  and  the  facts  that 

ln(ai  +  a2  +  a3)  >  i(ln  ai  +  In  a2  +  In  a3)  +  In  3 

O 

and  lim/<-_oo  £m[A0|T?]  =  1  ,  we  have 


K^ooj{ln^No^  ~  1  + 
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>  Vm*  l{^{lnpm+  ln[£m[iVol2]  -  1)}  +  ln[P.(2)]} 

m= 0 

>  iU(r7, 3) 


Then  by  substituting  from  the  upper  bounds  in  (20)  and  (22)  and  taking  the  limit,  we 
will  have 

Urn  iln{£?(tfo|*?]-  1  +  ^0?)} 

00  A 

<  lim  iln^JVoW-l  +  ^t?)} 

K—kx>  A 

The  reverse  inequality  also  holds  by  definition. 

Thus  the  same  threshold  can  be  used  for  all  sensors;  consequently  Ppt  =  Pf  and 
Pw.  =  PM  for  i-  —  1  Let  no  denote  the  number  of  times  that  the  fusion  center 

collects  the  decision  of  the  K  sensors  before  making  its  own  decision.  Then  in  analogy 
to  the  k  —  out  —  of  —  K  rule  of  block  data  fusion,  the  decision  test  of  the  fusion  center 
becomes 


decide  H\  if 
defer  " 
decide  Hq  " 


k(n0 )  >  6* 

6*  >  &(no)  >  a* 
a *  >  k(n0) 


where  fc(no)  is  the  number  of  sensors  which  have  transmitted  ‘1’  by  the  no-th  time  the 
fusion  center  has  collected  the  decisions  of  the  sensors  and  the  thresholds  ( a‘,b *)  are 
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induced  by  (a,  b )  and  are  functions  of  Pf,  Pm  and  K.  Then  Pf0  and  Pm0  have  forms 

Pf0  =  E  (fW(l -Pf)**-0  (25) 

<=f6-l  \  / 

Pa^o  =  F  E  (^)  ^(l  -  W  (26) 

Similar  to  (A2)  the  following  assumption  is  made  for  robustness  (see  Appendix). 
Assumption  4(A4):  For  sequential  data  fusion  Pp  <  \b*]/K  and  Pm  <  1  —  fa*]  j K  for 
all  K. 

Finally,  the  following  assumption  is  necessary  for  robustness  (see  [4]). 

Assumption  5(A5):  In  >■  Pm  In  z£e1, 

fmff 

Proposition  5:  Under  (Al),  (A3)  and  assumption  1  of  [2],  the  performance  index  defined 
by  (20)  satisfies 

*/<x>  (  Li ,  tjq,  t),  i  —  0,1,...) 

*^oo ( ^ *i  Vo,  Vi  i  —  1)  "  ') 

<  Joo{  Di,Vo,V,  *  =  o,  1,  — )  (27) 

Proof:  By  proposition  4,  the  same  threshold  can  be  used  by  all  sensors  to  obtain  asymp¬ 
totically  optimal  sequential  data  fusion.  Then  Pp0  and  Pm0  will  have  the  forms  described 
by  (22)  and  (23).  Therefore  a  similar  proof  to  that  of  proposition  2  can  be  developed 
corresponding  to  the  exponent  of  the  error  probability.  According  to  proposition  5  of  [4], 
the  least-favorable  pair  of  distributions  for  expected  sample  sizes  are  the  same  as  those 
for  the  error  probabilities  under  (A4).  Thus  the  saddle  point  of  performance  measure  (19) 
can  be  reached  by  the  same  least-favorable  pair  of  proposition  2  in  section  II. 
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IV.  CONCLUSION 


In  this  paper  we  formulated  and  solved  several  problems  of  robust  data  fusion  for 
situations  characterized  by  incomplete  knowledge  of  the  probability  distributions  of  the 
discrete-time  observations.  The  uncertainty  in  the  probability  distributions  was  modeled 
by  2-alternating  capacity  classes.  For  robust  block  data  fusion,  suboptimal,  optimal,  and 
asymptotically  optimal-as  the  number  of  sensors  increases-schemes  were  derived.  For 
sequential  data  fusion  an  asymptotically  robust  scheme  was  developed.  This  scheme  is 
robust  in  the  sense  that  the  sensors  use  likelihood  ratio  tests  and  the  fusion  uses  the 
sequential  probability  ratio  test  (SPRT).  In  all  situations  the  robust  fusion  rule  use  like¬ 
lihood  ratios  and  thresholds  that  depend  on  the  pair  of  least-favorable  distributions  in 
the  uncertainty  class.  The  advantage  of  asymptotic 'analysis  manifests  itself  in  that  the 
optimal  thresholds  of  the  various  sensors  turn  out  to  be  the  same,  which  in  turn  reduces 
considerably  the  required  computation  effort. 
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APPENDIX 


In  (A2),  it  follows  from  [2]  that 


f  =  [u,  g_  yg  -  PF)]-\m +Kln  i 

Pm  Pf  Pm 


and  for  large  A 


\k*}/K  «[ln^ — -%-K1  ~  £0  b  1  ~  Pf 

PmPf  Pm 

If  770  =  0,Pm  =  PF,  then  A  =  A/2  and  \k'}/K  *  1/2.  Thus  (A2)  is  not  restrictive. 


Similarly,  in  (A3) 


~  fy)]-'(°  +  A-  In  l^E 


PmPf 


Pm 


b-  =  Png--y(1  -  fF)]-»(t+  Jin. 


PmPf 


Pm 


and  for  large  A 


If  Pjvf  =  Pp,  as  A 


r«*l  »  [In  L1-  .  ^Af)(l  -  PF ),-i  ^  1  -Pf 
PmPf  Pm 

[A]  »  [In  0-~  Pm)(1  -  Pf),-i  ^  1-_Pf 
PmPf  Pm 

oo,  [a*] /A  «  1/2,  [A] /A  »  1/2. 


